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Abstract. In the framework of test-particle and cosmic-ray modified hydrodynamics, we calculate synchrotron 
emission radial profiles in young ejecta-dominated supernova remnants (SNRs) evolving in an ambient medium 
which is uniform in density and magnetic field. We find that, even without any magnetic field amplification 
by Raleigh- Taylor instabilities, the radio synchrotron emission peaks at the contact discontinuity because the 
magnetic field is compressed and is larger there than at the forward shock. The X-ray synchrotron emission 
sharply drops behind the forward shock as the highest energy electrons suffer severe radiative losses. 
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1. Introduction 

Shocks in supernova remnants (SNRs) are believed to pro- 
duce the majority of the Galactic cosmic-rays (CRs) at 
least up to the "knee" (~ 3 x 10 15 eV). The particle ac- 
celeration mechanism most likely responsible for this is 
known as diffusive shock accelera tion (DSA) (e.g., iDrurvl 
Il983t iBlandford fc Eichlerl Il987^ . This mechanism may 
transfer a large fraction of the ram kinetic energy (up 
to 50%) into relativistic particles and remove it from th e 
thermal plasma (see, for example. I Jones fc Eirisor]ri99llh 
Convincing observational support for the acceleration 
of particles in shell-type SNRs comes from their non- 
thermal radio and X-ray emissions due to synchrotron 
radiation from relativistic GeV and at least TeV elec- 
trons, respectively. In radio and X-rays, synchrotron- 
dominated SNRs display various morphologies: for in- 
stance, the synchrotron emission dominates in t wo bright 
limbs in SN 1006 (e.g., iRothenffug et alj[2004|) whereas 
it is distorted and complex in RX J1713. 7-3946 (e.g., 
ICassam-Chenai et alJl2004bj) . The detection and imaging 
with the HESS telescopes of TeV 7-rays in RX J1713.7- 
3946 provides unambiguous evidence for particle acceler- 
ation to very high energies. The 7-ray m orphology in this 
remn ant is similar to that seen in X-rays ( Aharonian et alJ 
120041) . 

Recent works based on Chandra llVink fc La ming 

l2003l for Cas A) and XMM-Newton 
l|Cassam-Chenai et alJ l2004al for Kepler's SNR) ob- 
servations have demonstrated that X-ray synchrotron 
emission is also present in ejecta-dominated SNRs and 
largely contributes to the continuum emission at the 



forward shock. This X-ray emission arises from sharp fila- 
ments encircling the SNR's outer boundary. The observed 
width of these filaments is a few arcseconds, and has been 
used to cons train the magnetic field intensity just b ehind 
the shock 1 l|Vink fc Laming] |2003t iBerezhko et alJ 12 003: 



. _ 13t L _ 

IBerezhko fc V61lJl2004 IVolk et al.ll2005l lBalledl2005T) 

A number of recent hydrodynamical models, includ- 
ing particle acceleration and photon emission, have been 
presente d to expl a in va rious features of these obser- 
vations. iRevnold s (1998) has described the morphol- 
ogy and spectrum of the synchrotron X-ray emission 
from SNRs in the Sedov evolutionary phase. Similar 
work based on numerical s i mulat ions was done by 
Ivan der Swaluw fc Achterberd l|2004f) who take into ac- 
count the diffusion of particles. CRs are treated as test- 
particles in these studies. 

Here, we expand on the work of IRevnolds! l|l998ji by 
considering young (ejecta-dominated) SNRs. We investi- 
gate the synchrotron emission morphology, both in radio 
and X-rays, as well as how it can be modified by efficient 
particle acceleration. Our results show that the radio and 
X-ray profiles are very different due to the effects of the 
magnetic field evolution and synchrotron losses in the in- 
teraction region between the contact discontinuity and the 
forward shock. For typical parameters, the radio emission 
peaks at the contact discontinuity while the X-ray emis- 
sion forms sheet-like structures at the forward shock. 



Magnetic field values are found to be at least 30 times 
higher than the typical Galactic field of 3/iG and imply that the 
field has been amplified, pe rhaps by the particle acceleration 
process llBell fc Lucekll200 J) . 
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2. Hydrodynamics and particle acceleration 

The hydro-dynamic evolution of young supernova rem- 
nants, including the backreaction from accelerated parti- 
cles, can be described by self-similar solutions if the initial 
density profiles in the ejected material (ejecta) and in the 
ambient m edium have power-law distributions l|Chevalierl 
Il982lll983ll and if the acceleration efficiency (i.e. the frac- 
tion of total ram kinetic energy going into suprathermal 
particles) is independent of time. 

Here, we use the self-similar model of lChevalierl l)l983l) 
which considers a thermal gas (7 = 5/3) and the 
cosmic-ray fluid (7 = 4/3), with the boundary condi- 
tions calculated from the non-linear diffusive s hock a ccel- 
eration (PSA ) mod el of iBerezhko fc Ellison] dl999l) and 
Elli son et alJ l|2000h as described in iDecourchelle et alJ 
(2000). This acceleration model is an approximate, semi- 
analytical model that determines the shock modification 
and particle spectrum from thermal to relativistic ener- 
gies in the plane-wave, steady state approximation as a 
function of an arbitrary injection parameter, ry; n j (i.e. the 
fraction of total particles which end up with suprather- 
mal energies) . The validity of the self-similar s olutions has 
been discussed by IDecourchelle et alJ l|2000|) and direct 
comparisons between this sel f-similar mod e l and the more 
general CR- hydro model of lEllison et alJ l|2004|) showed 
good correspondence for a range of input conditions. 

The hydrodynamic evolution provides the shock char- 
acteristics necessary to calculate the particle spectrum at 
the forward shock 2 , at any time. Once a particle spectrum 
has been produced at the shock, it will evolve downstream 
because of radiative and adiabatic expansion losses. We as- 
sume that the accelerated particles remain confined to the 
fluid element in which they were produced, so adiabatic 
losses are determined directly from the fluid element ex- 
pansion. The basic power law spectrum produced by DSA, 
before losses are taken into account, is modified at the 
highest energies with a exponential cutoff, exp(— p/p max ), 
where p max is determined by matching either the accel- 
eration time to the shock age or the upstream diffusive 
length to some fraction of the shock radius. In our simu- 
lation, the electron-to-proton density ratio at relativistic 
energ ies, (e/p) Te \, is set equal to 0.01 (see lEllison et all 
2000). 

Unless explicitly stated, our numerical examples are 
given for the following supernova parameters: M e j = 5 M 
for the ejected mass, -E51 = 1 where -E51 is the kinetic 
energy of the ejecta in units of 10 51 ergs and n = 9, where 
n is the index of the initial power-law density profile in 
the ejecta (p cx r~ n ). In our simulations, the SNR age is 
tf = 400 years and the shock velocity at the forward shock 
is v s ~ 5 x 10 3 km/s. For the ambient medium parameters, 
we take a magnetic field Bq = 10 pG, a density no = 

2 We do not consider CR production at the reverse shock 
since the magnetic field at the reverse shock may be consider- 
ably smaller than that at the forward shock due to the dilution 
by e xpansion and flux f reezing of the progenitor magnetic field 
(see lEllison et af]l2005l) . 



0.1 cm -3 , an ambient gas pressure p & fl/k = 2 300 Kcm -3 
and s = 0, where s is the index of its initial power-law 
density profile (p oc r~ s ). The case s = corresponds to 
a uniform interstellar medium (s = 2 describes a stellar 
wind). 

In the next section, we discuss the importance of the 
magnetic field for the synchrotron emission and particle 
acceleration. We do not, however, explicitly include the 
dynamical influence of the magnetic field on the hydrody- 
namics. 

3. Results 

3.1. Magnetic field 

To track the synchrotron losses, we are interested in the 
temporal evolution of the magnetic field behind the shock. 
We assume the magnetic field to be simply compressed at 
the shock and passively carried by the flow, frozen in the 
plasma, so that it evolves conserving flux. In this sim- 
ple 1-D approach, we do not consider any production of 
the SNR magnetic field, for instance, by hydrodynami- 
cal instabilities which is an additional effect. As for the 
magnetic field ahead of the forward shock, it is assumed 
to be isotropic and fully turbulent. Appendix E] (see the 
on-line version) shows how to compute the magnetic field 
profile for self-similar solutions in both test-particle and 
nonlinear particle acceleration cases. 

3.1.1. Test-Particle limit 

We first discuss the behavior of the normal and tangen- 
tial components of the magnetic field in the test-particle 
case where the backreaction of the accelerated particles is 
neglected. 

When the SNR evolves in an ambient medium which is 
uniform in density and magnetic field, the expansion and 
flux freezing generally cause the tangential component of 
the magnetic field to increase at the contact discontinuity 
whereas the normal component falls to zero (Fig. ^) . As 
a result, the magnetic field profile is dominated by the 
tangential component. 

One has often invoked hydrodynamic instabilities 
to explain the magnetic field increase at the interface 
between the shocked eje cta and the shocked ambient 
medium ll.Iun et all Il99.il The numerical simulations of 
ljun &: Normanl l(l996|) have shown that the magnetic field 
could be amplified by a factor 60 by Rayleigh- Taylor and 
Kelvin-Helmholtz instabilities. Here, we note that simple 
advection of the magnetic field already predicts amplifica- 
tion by a factor 5 (Table |A~T1 top, n = 9). 

We note that, if the SNR evolves in a wind with a de- 
creasing initial density profile, advection goes the other 
way (diluting the magnetic field instead of amplifying it). 
But when both the ambient density and magnetic field 
decrease with radius, as would be the case for a pre- 
supernova stellar wind, the magnetic field is larger close 
to the contact discontinuity than at the forward shock 
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Fig. 1. Radial profile of the normal (B r ) and tangential 
(B t ) components of the magnetic field in a test-particle 
self-similar model. Each component is normalized to the 
forward shock. 
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Fig. 2. Magnetic field radial profile for different values of 
the injection efficiency, ryi n j, when the hydrodynamics is 
coupled with the non-linear DSA model. The width of the 
shocked region is smaller and smaller as the feedback of 
the accelerated particles on the SNR dynamics increases. 
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Table 1. Compression ratio, r to t, and downstream mag- 
netic field, B s , at the forward shock obtained for differ- 
ent injection, ?7i n j. The magnetic field compression ratio is 
given by r B = B s /B a = ^1/3 + 2 r t 2 ot /3 and B = 10 M G 
here. 
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Fig. 3. Radio (top panel) and X-ray (bottom panel) syn- 
chrotron volume emissivity, e„, radial profile for different 
injection efficiencies. 



(by a factor of ~ 1000 in some cases). This is because 
the dilution of the advected magnetic field is negligible 
compared to the fact that the ambient magnetic field was 
much larger at early times. 

3.1.2. Nonlinear Particle Acceleration 

We now consider the behavior of the normal and tangen- 
tial components of the magnetic field in the nonlinear case 
where the backreaction of the accelerated particles on the 
shock is taken into account. 

In the ideal non-linear case, where the acceleration is 
instantaneous, the magnetic field diverges at the contact 
discontinuity because of its tangential component, what- 



ever the injection efficiency is, as in the test-particle case. 
However, the contrast between the magnetic field in a 
given fluid element and the one just behind the shock, 
will be always smaller than in the test-particle case (see 
Table IA.1(I . Figure shows the profile of the total down- 
stream magnetic field for different values of the injection 
efficiency. Tabled shows the associated compression ratio 
and immediate post-shock magnetic field. 

3.2. Synchrotron emission 

Once the magnetic field structure and the particle spec- 
trum (attached to a fluid element) modified by the ra- 
diative and adiabatic expansion losses as computed in 
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onto the line-of-sight for different injection efficiencies. 
Emission from particles accelerated at the reverse shock 
is assumed to be negligible. 



iRevnoldsl lll998l) are known, we comp ute the synchrotron 
emission ijRvbicki k, Lightmanlll979|) . averaged over the 
pitch-angle, in any energy band 3 . 

Figure shows the radial profiles of the synchrotron 
emission in the radio (top panel) and X-ray (bottom 
panel) domains for different injection efficiencies, r^j. An 
increase in the injection efficiency not only provides a 
larger number of accelerated electrons, but also a larger 
compression of the downstream magnetic field (see Table 
^) and a narrower interaction region. These effects com- 
bine to produce enhanced synchrotron emission as the in- 
jection increases. 

The radio synchrotron emission is produced by GeV 
electrons which are not affected by radiative losses. 
Consequently, the radio synchrotron emission critically 
depends on the final magnetic field profile (Fig. EJ) and, 
therefore, peaks at the contact discontinuity. In contrast, 
the X-ray synchrotron emission is produced by the highest 
momentum electrons (~ 10 3-5 m p c) which, depending on 
the downstream field strength, may suffer radiative losses. 



We did not calculate the synchrotron emission from the 
precursor. 



The high energy electrons that have been accelerated at 
the earliest time have suffered strong synchrotron losses 
as they were advected behind the shock. Because of this, 
they are not numerous enough at the end to radiate in 
the X-ray regime despite a strong magnetic field. As a 
result, the X-ray synchrotron emission rapidly decreases 
behind the shock. The X-ray profile becomes sharper when 
the injection efficiency increases because it provides larger 
compression of the downstream magnetic field and then 
stronger synchrotron losses. 

Figure 0] shows the synchrotron emission after inte- 
gration along the line-of-sight. The radial profile of the 
radio emission (top panel) shows a peak at the contact 
discontinuity. The radial profile of the X-ray projected 
synchrotron emission (bottom panel) shows bright rims 
just behind the forward shock whose width decreases as 
the injection efficiency increases. 

4. Discussion and Conclusion 

We have computed the radio and X-ray synchrotron emis- 
sion in young ejecta-dominated SNRs. This has been done 
using a one dimensional, self-similar hydrodynamical cal- 
culation coupled with a non-linear diffusive shock accel- 
eration model, and taking into account the adiabatic and 
radiative losses of the electron spectrum during its advec- 
tion in the remnant. 

We show that the morphology of the synchrotron emis- 
sion in young ejecta-dominated SNRs is very different 
in radio and X-ray. This is the result of the increased 
magnetic field toward the contact discontinuity, to which 
only low energy electrons that emit radio are sensitive, 
while the high energy electrons emitting X-rays experi- 
ence strong radiative losses and are mostly dependent on 
the post-shock magnetic field. 

Briefly, the radio synchrotron emission increases as one 
moves from the forward shock toward the contact discon- 
tinuity due to a compression of the magnetic field (par- 
ticularly its tangential component), assuming both uni- 
form ambient density and upstream magnetic field. Such 
a compression naturally results from the dynamical evo- 
lution of the SNR. In contrast, because of the radiative 
losses, the X-ray synchrotron emission decreases behind 
the forward shock and forms sheet-like structures after 
line-of-sight projection. Their widths decrease as the ac- 
celeration becomes more efficient. 

The morphology of the radio synchrotron emission ob- 
tained for the young ejecta-dominated stage of SNRs will 
differ from tha t of SNRs in the Sedov phase (but not in 
X-ray). Indeed. IRevnoldsl l|l99ct) has shown that both the 
normal and tangential components of the magnetic field 
decrease behind the forward shock in the Sedov phase and, 
as a result, we expect the radio synchrotron emission to 
decrease behind the shock (however, less rapidly than the 
X-ray synchrotron emission since the radio electrons do 
not experience radiative losses). 

Our model qualitatively reproduces the main features 
of the radio and X-ray observations of emission in young 
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ejecta-dominated SNRs (e.g., Tycho and Kepler), i.e. 
bright radio synchrotron emission at the interface between 
the shocked ejecta and ambient medium, and a narrow fil- 
ament of X-ray emission at the forward shock. However, 
this model is unable to reproduce the thin radio filaments 
observed at the forward sho ck in some SNRs ( for instance 
those seen in Tycho's SNR. iDickel et al.lll991l) . 

We note that extensions of this work to cases with 
exponential ejecta profiles and/or SNRs evolving in a pre- 
supernova stellar wind with varying magnetic fields, can- 
not be done with self-similar solutions. These cases can be 
calculated in the n umerical CR-m odificd hydrodynamical 
model descr ibed in lEllis on et al J 1120051) and this work is 
m progress l|Ellison fe Cassam-Chenaill2005[) . 
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Appendix A: Magnetic field evolution 

The evolution of the normal (subscript r) and tan- 
gential (subscript t) components of the magnetic 
field at the downstream position, B, is given by 
( Reynolds fc Chevaherlll98ijk 



B r (r) 



Bt{r) = Bt. 



Br,j — 

p r 



Pj 



(A.l) 
(A.2) 



and the total magnetic field is simply l|Revnoldslll99^ : 



B(r) = (B r (rf + B t (r) 2 ) 



2\!/2 



(A.3) 



In these equations, r and p are, respectively, the radius 
and density of a fluid element at the current time that 
was shocked at the previous time tj . At time tj , the fluid 
element was just behind the shock at the radius Tj , with 
a density pj and a magnetic field Bj . 

We assume that the upstream magnetic field at time 
tj, Bo j is isotropic and fully turbulent so that the 
components of the immediate post-shock magnetic field 
Bj in Eas HA.1|) an d ljA.21) are given on average by 
l|Berezhko et al.ll2002() : 



B r j — 1/\/SBq 



B. 



/ 2/3 r to t Bqj 



(A.4) 
(A.5) 



where r tot is the shock compression ratio. In the self- 
similar approach, rw is assumed independent of time (see 
iDecourchelle et ai1l2000t for details). 

We consider that the current magnetic field upstream 
of the forward shock, Bq, s , can behave like: 



O.s 



Bo. 



(A.6) 



where r s is the current shock radius. If the magnetic field is 
uniform, the index q is equal to 0. In a stellar wind (s = 2), 
the magnetic field profil e may be decreasing yielding q = 1 
I Lvutikov fc Pohill2004h or q = 2 if we assume that it is 
frozen in the plasma. 

We define the magnetic field contrast factor, <tb = 
B/B s , as the ratio between the current magnetic field in 
a fluid element, B, and the current one just behind the 
shock, B s . We have: 



2r t 2 ot 



a 



B, 



1 



2r t 2 ot 



1/2 



(A.7) 



where ob t = B r /B r s and ag t = B t jB t s are the magnetic 
field contrast factors of the normal and tangential compo- 
nents of the field, respectively. The components B r s and 
B t . s obey the same relation as in Eqs l|A.4|) and (|A.5|I . 



A.l. Test-Particle limit 

Assuming adiabaticity of the thermal gas, the magnetic 
field contrast factors of the normal and tangential compo- 
nents of the field are given by: 



<JB r = 



Rs 
R 
p 

P. 



-3/5 



0r 

R 



-(ll-3s)/5 



where the indexes (3 r and fit are given by: 

- 3 



r = (q-2) 



(it = 



3-s 
33 — 'is(n — 5) 



5(3 -s) 



n — 3 
3- s' 



(A.8) 
(A.9) 

(A.10) 
(A.ll) 



In Eqs HA.8|) and (|A.9|) . R s /R and P g , s /P g are the ratio 
of the self-similar radii and thermal gas pressures, respec- 
tivel y, between the s hock (subscript s) and a fluid element 
Csee Ichevalierlll982|h They depend on n and s, but also 
weakly on Vj/v s where v s and Vj are the current shock ve- 
locity and the shock velocity at the time tj, respectively. 

In the framework of these self-similar solutions, the 
forward shock velocity tends to infinity at early times, cor- 
responding to fluid elements close to the contact disconti- 
nuity at the current time. To limit the maximum velocity 
to a realistic value, we look at the value of gb for a shock 
velocity ratio Vj/v s = 10. For the typical forward shock 
velocity v s that we have used for the numerical applica- 
tion, the initial velocity corresponds to d,- ~ 5 X 10 4 km/s. 
This shock velocity is the criterion used to define the ra- 
dial position of the oldest fluid element that is currently 
located close to the contact discontinuity. 

Here, we consider the case of both an uniform ambient 
medium (s — 0) and upstream magnetic field (q = 0). 
Under this assumption, B s = Bj, since r to t is constant 
with time. Then, the magnetic field contrast factor, <jb is 
equal to B/Bj and can be viewed as a compression or a 
dilution factor. Table fA~TI ftop) gives the contrast <jb for 
different values of n. 



A.2. Nonlinear Particle Acceleration 

In the ideal non-linear case, where the acceleration is in- 
stantaneous and efficient, the thermal gas pressure falls to 
zero at the contact discontinuity while the relativistic gas 
pressure goes to infinity. Hence, the contrast factor of the 
tangential field component, <jB t , given by Eq. (|A.9|I . ob- 
tained in the test-particle limit, is not defined when Vj/v s 
tends to infinity. 

However, the contrast of the tangential component of 
the magnetic field can also be found by using the adia- 
baticity of the relativistic gas: 



P 



C,5 



-3/4 



-(10-3s)/4 



0t 



(A.12) 
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n 


Rs/R 






OB 


v a (km/s) 




7 


1.181 


0.964 




0.95 


4840 


Test-Particle 


9 


1.140 


0.885 




5.0 


4850 




12 


1.121 


0.836 




54 


4940 




n 


Rs/R 


P s /P P c , s /P 


min(o\B 


) max(CTs) 


v s (km/s) 




7 


1.080 


1.135 0.806 


0.50 


0.66 


4370 


Nonlinear DSA 


9 


1.060 


1.045 0.754 


2.6 


3.4 


4470 




12 


1.051 


0.988 0.714 


28 


36 


4610 



Table A.l. Test-Particle: Magnetic f ield contrast fac tor, erg, for a velocity ratio Vj/v s — 10, radius ratio and self- 
similar thermal gas pressure ratio (see IChevalieiJll982h and shock velocity, for different values of the indexes n (s = 
and q = 0). Nonlinear DSA: same as Test-Particle but with lower and upper limits on the magnetic field contrast 
factor, <7b, with i] in j = 10 -3 (see Sect. IA.2l for more explanations). 



where the index f3' t is given by: 

4n-30-3 S (n-5) +g n-3 
4(3 — s) 3 — s 

In Eq. (|A.12(1 . P c . s /P c is the ratio of the self-similar rela- 
tivistic gas pressures between the shock (subscript s) and 
a fluid clement. This ratio depends on n, s, and Vj/v s . 
The contrast of the normal field component, <7B r , is still 
given by Eq. ljA.81) . The asymptotic behavior of the con- 
trast factor, (7 q ± , can be derived from Eq. (|A.12|I because 
the relativistic gas pressure does not tend to zero at the 
contact discontinuity. 

Because the thermal gas pressure vanishes as we ap- 
proach the contact discontinuity in the case of ideal par- 
ticle acceleration, i.e., when the acceleration is instanta- 
neous and efficient, the contrast of the tangential field 
component, as t , will always be smaller than in the test- 
particle case where the thermal gas pressure rapidly tends 
to a constant (see Eq. IA.9|1 . Table IA.1I (bottom) gives 
the lower and upper limits on the magnetic field con- 
trast factor, o\b, in the case of ideal nonlinear particle 
acceleration for ?7i n j = 10~ 3 and for different values of 
n when both the ambient medium and upstream mag- 
netic field are uniform (s — and q — 0). The lower 
and upper limits on ub are obtained by replacing in Eq. 
(|A. 12|) the ratio of the self-similar relativistic gas pres- 
sures, P c . s /P c , by the ratio of the self-similar total gas 
pressures, P s /P = (P c , s + P S , S )/{P C + P s ), an d by the 
ratio between the self-similar relativistic gas pressure at 
the shock and the self-similar total gas pressure, P c . s jP, 
respectively. 

However, for an injection efficiency lower than ~ 5 x 
10~ 4 , the acceleration is not efficient enough for the shock 
to be modified at the beginning of the evolution. In that 
case, the fluid elements that have been shocked at the 
earliest times are still dominated by the thermal gas so 
that test-particle solutions could still apply locally. 



